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The theory of the strong interactions, Quantum Chromodynamics (QCD), has been 
addressed by a variety of non-perturbative techniques over the decades since its 
introduction. We have investigated Hamiltonian formulations with different quan- 
tization methods and approximation schemes. In one method, we utilize light-front 
coordinates to investigate the role of bosonic zero modes in leading to confinement. 
In another method we are able to obtain spectra for the mesons and baryons using 
constituent quark masses but no phenomenological confinement. We survey our 
principal accomplishments to date and indicate our future directions. 



1 Introduction and Motivation 

For strong interaction physics, one might well ask 'Why develop Hamiltonian 
methods for gauge theories - after all, does not the lattice gauge method work 
well?' 

While lattice gauge methods work well for certain observables, we have 
a wider range of observables in mind. In addition, we are motivated by the 
desire to work within a Hamiltonian framework which we find more physically 
intuitive. One final motivation for us is that we believe there are potential 
advantages of using advanced methods from quantum many-bdoy theory. 

On the other hand, we are also quick to acknowledge many challenges 
which we face. For example, little is known about renormalization and scale 
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dependence within the Hamiltonian approach to strong interactions. We will 
lose manifest gauge invariance when we derive a Hamiltonian expressed only 
in terms of the independent degrees of freedom. Any approximations, such 
as a truncation, will usually lead to gauge-dependent results. Finally, due to 
our approximations, we may lose other symmetries respected by the original 
Lagrangian. 

Our purpose here is to outline our recent progress in implementing Hamil- 
tonian approaches and to indicate where we still face major hurdles. 

2 Introduction to Many-Body Theory of H e ff 

Here we discuss certain aspects of recent developments in the theory of effective 
Hamiltonians for quantum many-body systems which are particularly relevant 
for our applications to strong interactions. Note that these developments are 
cast in a form which is independent of the kinematics of the interacting parti- 
cles. 

The basic framework we adopt has been utilized extensively in nuclear 



for the eigenenergies E a and the eigenstates |* a ) of the many-particle system, 
where a is some label characterizing the states. But it is impossible to solve this 
problem in the full Hilbert space S when the number of particles in the system 
exceeds 3 or 4 because it contains too many degrees of freedom. Consequently, 
one wishes to truncate the problem to a smaller space S of dimension d, in 
which it becomes tractable to carry out the calculation. Now let |$^) represent 
the projections of d of the states into S. Thus we define the effective 

Hamiltonian H e f f in S to satisfy 



where the eigenvalues {Ep} are d of the exact eigenvalues {E a } in Eq. (|lj) . Be- 
cause the are projections of the \^ a ), they are, in general, not orthogonal. 

The question then arises whether an appropriate H e f / exists for any given 
truncation. One can show this to be true by constructing the biorthogonals 
to l^), namely, |$ 7 ), which satisfy ($ 7 |<i> / 3) = <5 7i a. It then follows that the 
effective Hamiltonian H e ff always exists and is of the form 




H eff \$ l3 ) = E0\$0), 



(2) 




(3) 



pes 
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which automatically satisfies Eq.(|^). As KirsonB has emphasized, the question 
is not whether H e ff exists, but whether it has a simple enough form, so as to 
be useful. 

We use the time-independent-perturbation-theory approach!-!!!'! in estab- 
lishing the connection between -ff e // and H. The basic idea involves the sepa- 
ration of the Hilbert space into two parts, using the projection operators P and 
Q, where P defines the truncated or 'model' space, defined by the eigenstates 
of an unperturbed Hamiltonian Hq, and Q defines the excluded space outside 
the model space. The projection operators P and Q define non-overlapping 
spaces, so that PH Q = 0. 

In the full Hilbert space, a typical many-body choice for H is of the form 

A A 

H = ^2 ti + 51 Vij = T + V = (T + U) + (V — U) = H + Hi, (4) 

i—l i<j 

where U is some single-particle or 'auxiliary' potential, Hq = T + U , and 
Hi = V — U is the residual interaction. Only two-body interactions have 
been assumed among the A-particles, but the method can be generalized to 
many-body forces. _ 

Using the Feshbach projection methodQ, one can explicitly project H into 
the P and Q spaces and rewrite the P space equation (omitting the subscript 
a everywhere) in the form 



PHP + PHQ E^QHQ QHP 



P\*)=EP\V), (5) 



where P\^) = \$). The term in the square brackets defines the effective 
Hamiltonian 

H eff =PH P + V{E), (6) 

where 

V(E) = PHiP + PHiQ - _ l QHQ QHiP (7) 

is the effective interaction. It should be noted that, in general, V is an A- 
particle operator and the energy E in the denominator corresponds to one of 
the exact eigenenergies of the A-particle system. 

Although V(E) is an ^4-particle interaction, the standard assumption is 
to approximate it in terms of a perturbation-theory expansion using two-body 
interactions. There are many uncertainties-associated with this approach for 
conventional nuclear physics applications t3 and experience gained in those 
applications may well assist us in our applications to quantum field theories. 
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Our approach is to take the ^4-particles as all active (i.e. we do not assume 
a passive 'core'), and Eq.(|^) may be interpreted as a generalized ^4-particle G- 
matrix equation. For a one-dimensional model space, the exact solutions for 
the eigenvalues are given by 

E = E Q + V(E), (8) 

where Eq is the eigenenergy of the unperturbed Hamiltonian PHqP. If V(E) 
can be constructed for the A-aaxticle system, then Eq.(||) can be solved di- 
agrammatically or iteratively em, to obtain all the eigenenergies of the A- 
particle system whose eigenstates have non- vanishing projections on the one- 
dimensional model space. Examples of these procedures with soluble models 
prove instructive Ej. 

It is not generally possible to construct the fulLd-particle G matrix. We 
approximate it by the two-particle reaction matrix Gliil plus higher-order terms. 
The two-particle G matrix is simply the infinite sum of two-particle ladder 
terms. 

The perturbation-theory expansion for V{E) is now rewritten as a per- 
turbation series in G[uS) where u> is an arbitrary energy, the 'starting' energy, 
around which we make an expansion. In our application to no-core model 
spaces, these corrections at the two-particle level are all of the folded-diagram 
type. The remaining corrections generate effective many-body interactions 
which are neglected in the applications we discuss here. 



3 Renormalization of H e ff 

In this section, we introduce a way of implementing WilsoKurenormalizationlll 
within the context of the theory of effective Hamiltonians.E-2 Our renormaliza- 
tion scheme involves manipulations at the level of the generalized G-matrix. 
We show how to calculate the beta function within this context and exhibit 
our method using simple scale-invariant quantum mechanical systems. 

We have argued above that the knowledge of the matrix G allows us to 
obtain H e ff. In what follows, we shall therefore restrict our attention only to 
G. For the sake of convenience we choose to work in the momentum repre- 
sentation where the kinetic energy term in the Hamiltonian is diagonal. To 
introduce the concept of renormalization we shall focus our attention on the 
one-particle system. The matrix elements of G are here given by 

G kk , = (k\PVP\k>) + 

+ f dpdp'{k\PVQ\ P ){p\ ^—\ p ')(p'\QVP\k')+--- (9) 
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where we have set U = and expanded QVQ out of the denominator. 

Let us suppose that the potential V depends on a single coupling constant 
Ho, which we shall call the bare coupling constant. It is clear from Eq. (||) 
that the matrix element Gkk' will be a function of fio- The expression in 
Eq. @ may, in general, require regularization due to the divergence arising 
from the integral. The regularization that we choose consists of introducing 
an ultraviolet cutoff A. The matrix element in Eq. (||) is now a function of the 
coupling constant fio and the cutoff A. At the end of the calculation we must 
remove the cutoff, i.e. we must take A to oo, which, as discussed above, may 
in general lead to divergence. One way to avoid the divergence is to replace 
the coupling constant fj-o with a function of A, which we denote as /x(A), and 
then require that matrix element in Eq. (^) remain finite and independent of 
the cutoff as the cutoff is removed. In other words, we demand that 

lim -^GW(A,/i(A)) = (10) 

The function //(A) thus plays the role of the renormalized coupling constant. 

The dependence of the coupling constant on the cutoff is usually expressed 
in terms of the beta function, which is defined by 

«rt-A*. (ID 

Within our formalism, Eqs. (|ic|) and (O) can be used to calculate the beta 
function. 

Note that once Eq. ( fL0|) is satisfied and /u(A) is determined, then H e g, 
based on Gkk' (A, n(A)), should also be independent of A as A — > oo. Thus, 
the complete problem of renormalization is solved. 

We shall now illustrate the method prescribed above in two simple cases of 
a Dirac particle in 1 dimension and a Schrodinger particle in 2 dimensions 113113. 
In both these cases the interaction potential will be taken as a delta function 
in position space : 

V{x) = -^ n \x), (12) 

where n is the dimension of configuration space. In the momentum space the 
interaction potential would simply be a constant, i.e., 

V(k) = -no. (13) 

We will choose Hq to be the pure kinetic operator, and our model space to 
consist of all states with momenta less than A. Thus Q projects onto the 
momentum range [A, oo]. 
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With the choice of the interaction potential described above, the series in 
Eq. (0) can be summed exactly and is given by 



G kk > = — % -M k-k% (14) 

1 + p I(UJ) 



where is given by 



J(w) ee / <Tp — -, (15) 

Jx uJ-Eo(p) 

Following the preceeding discussion we now introduce an ultraviolet cutoff A. 
Replacing po by the renormalized coupling constant p and using Eqs. ( |l0| ) and 
(|i"i"l), we obtain the beta function as 

8T 

m = M 2 A— . (i6) 

To obtain the explicit expression for the beta function we need to evaluate 
the integral appearing in Eq. (|l5|). For the 1 dimensional Dirac particle we 
have n = 1, Eq(p) — p + m and 

= / dp — — — = - In ) . (17) 

The corresponding beta function is given by 

f3(p) = -p 2 . (18) 

For the Schrodinger particle in 2 dimensions we have n — 2 and Eo(p) = 
p 2 /2 (we set the mass of the particle to unity.) Proceeding exactly as before, 
we obtain 

J( W ) = -27rln (^p) (19) 

and 

(3 = -Anp 1 (20) 

Note that the results in both examples above have the desirable property 
that the beta function is independent of the model space cutoff, A. The beta 
functions calculated give rise to asymptotically free theories and generate the 
generally accepted pattern for the flow of the coupling constant for the two 
examples described above. 
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4 Light-Front QCD Application 

Rather than take a pure light-front quantization approach, we opt for quantiz- 
ing on a space-like surface in 'near light-front coordinates' in order to maintain 
contact with many of the known results from the usual equal-time quantization 
method. B 

Here we will sketch the dynamics of the gluonic zero modes of the color 
SU(2) QCD Hamiltonian and show how the strong coupling solutions serve 
as a basis for the complete problem. Only an external charge density p m is 
considered here. 

4-1 Color SU (2) QCD Hamiltonian 

The Lagrangian in the near light front coordinate system reads 

C = ln_F%_ + £ (f^FI, + ^F^FU) - \f&F& - p a m A a + , (21) 

where the color index a is summed from 1 to 3, and the transverse coordinates 
are labeled by i — 1,2. We will also use the matrix notation; for example 
A- = A"t°/2. The field strength tensor contains the commutator of the 
gauge fields and the coupling constant g: 

F^v = d^A v - d v A^ - ig[A^ A v ] . (22) 

The A+ coordinates have no momenta conjugate to them. As a con- 
sequence, the Weyl gauge A" = is the most natural starting point for a 
canonical formulation. The canonical momenta of the dynamical fields At , A^ 
are given by 

F$_, 

F-i + W 2 F+i ■ (23) 
^From this, we get the Weyl gauge Hamiltonian density 

n w = \w_ w + \f&f& + -L J2 ( n ? - F -i) 2 ■ ( 24 ) 

Z Z LT] i=l,2 

We choose periodic boundary conditions in x~ = x ^ and xj_ on intervals 
of size [0,L]. Using the appropriate periodic delta functions, the quantization 
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n? = 



dC 
dFl_ 

dC 
9FL 



is straightforward. However, the Hamiltonian has to be supplemented by the 
original Euler-Lagrange equation for A + as constraints on the physical states 

G a (x ± ,x-)\<f>) = (D a _ b n b _ + Df ni + gp a m ) |$) 

= (D ab lt_ + G\) |$> = 0, (25) 

with the covariant derivatives 

D at _ g_s ab + gf acb A c _, 

Df = d ± S ab + gf acb A c L , (26) 

where f acb are the structure constants of SU(2). These equations are known 
as Gauss' Law constraints. Since the Gauss' Law operator commutes with the 
Hamiltonian 

[G a (x±,x-),H w ]=0, (27) 

time evolution leaves the system in the space of physical states. Furthermore, 
H\y is invariant under time independent resicLual gauge transformations whose 
generator is closely connected to Gauss' Lawlia. 

In order to obtain a Hamiltonian formulated in terms of unconstrained 
variables, we eliminate the A~ . Classically this would correspond to the light 
front gauge A- = 0. However, this choice is not compatible with gauge in- 
variance and periodic boundary conditions. Only the (classical) Coulomb light 
front gauge (<9_A_ = 0) is legitimate. The reason is that A- carries infor- 
mation on gauge invariant quantities, such as the eigenvalues of the spatial 
Polyakov (Wilson) loop 



V(x±) — Pexp 



ig / dx A_ (x± , x ) 



(28) 



which can be written in terms of a diagonal matrix a_(xjj 

)] V* . (29) 

Thus, we obviously need to keep these 'zero modes' a_(xj_) as dynamical 
variables. In order to eliminate the conjugate momentum, n_ of A—, by 
means of Gauss' Law, one needs to 'invert' the covariant derivative After 
an unitary transformation _D_ simplifies significantly (compare to Eq. (|26|)) 

L>_ =cL -ig[a-, . (30) 

Now Gauss' Law can be readily resolved: in the physical space one can make 
the replacement 

U-(x±,x-) ^p-(x±)- (dZ 1 )G ± {x ±7 y-). (31) 



8 



The zero mode operator p_(x±) is the conjugate momentum to a,-(x±). It 
has eigenvalue zero with respect to gL, i.e. d-P- = 0, and is therefore not 
constrained. 

The appearance of the zero modes also implies residual Gauss' Law con- 
straints. In the space of transformed physical states |x), they read 

J dx-Gi \ X ) = J dx- (Dfu'i + gp 3 m ) \ X ) = 0. (32) 

These two-dimensional constraints can be handled in full analogy to QED, 
since they correspond to the diagonal part of color space. This further gauge 
fixing in the SU(2) 3-direction is done via another gauge fixing transformation, 
which leads to the Coulomb gauge representation in the transverse plane for the 
neutral fields. In other words, we eliminate the color neutral, x~ -independent, 
two-dimensional longitudinal gauge fields 

a l ±(x ± ) = \l dy-dy^x^ - y ± )V± (V± ■ A 3 ± (y ± ,y_)) T - . (33) 

Here we use the periodic Greens function of the two dimensional Laplace op- 
erator 

L n^6 Pn L 

where ft = {n\,n2) and n\,n<i are integers. The conjugate momenta of these 
fields, p i L (x±), are defined analogously. Resolution of the residual Gauss' Law 
allows one to replace them, in the sector of the transformed physical space 
1$'), by the neutral chromo-electric field 



e ± (x ± )=gV ± j dy-dy ± d(x ± -y ± ){f 3ab A a L (y ± ,y-)U b L (y ± ,y-) 

T 3 

+p 3 m (y±,y~)}Y ■ ( 35 ) 

It is convenient to introduce the unconstrained gauge fields and their conjugate 
momenta: 

A / _ L (x ±7 x-) = Aj_(x ± ,x-) - a,j_(x±), 

W ± (x±,x~) = Ilj_(xj_,x-) - -p e ± (xj_). (36) 

These relations turn out to be important for neutral gluon exchange; recall 
that the subtracted fields are diagonal in color space. Note that the physical 
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degrees of freedom A' ± and W ± still contain (x±, x~)-indepcndcnt, color neu- 
tral, modes. Therefore, there is a remnant of the local Gauss' Law constraints 
- the global condition 



Its physical meaning is that the neutral component of the total color charge, 
including external matter as well as gluonic contributions, must vanish in the 
sector of physical states. 

The final Hamiltonian density in the physical sector explicitly reads 



where p and q are matrix labels for rows and columns, a_ g = (a_) gg and 
the prime indicates that the summation is restricted to n ^ if p = q. The 
operator G' ± (x± ,x~) is defined as 



4-2 Zero Mode Dynamics 

The principal advantage of an exact light front formulation is the apparent 
triviality of the ground state which simplifies calculations of the hadron spec- 
trum. The light front vacuum, however, is not guaranteed to be trivial in the 
zero mode sector. 

As can be seen from the dispersion relation for massless particles on the 

2 

light front, p + = soft momentum modes become high energy states. In 
this way, high energy physics becomes tied to long range physics, contrary to 
the equal time formulation. This physics appears in deep inelastic scattering 
at small scaling variable and is related to the long distance features of the 
proton. We will focus on the zero mode sector in order to try to acquire some 
insight into its dynamics. 

^From the comparison of abelian and non-abelian theories, striking differ- 
ences show up in the zero mode sector. Recently, in the equal time formalism, 




(37) 





(39) 
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the zero mode sector in QCD has been claimed to be relevant for the confine- 
ment phenomenon!^. On the level of approximations and restrictions followed 
below, the formal differences between light front and equal time approach are 
rather small and, consequently, results and methods are similar. 

In this work we do not restrict ourselves to the strong coupling approx- 
imation. We will, however, start with the strongly coupled theory to define 
our set of basis functions. We will restrict ourselves to pure gluonic SU(2) 
Yang-Mills theory. It already has the typical non-abelian features such as the 
Coulomb term which explicitly contains the zero modes in the denominator 
and the non-standard kinetic energy for the zero modes. 

The zero mode degrees of freedom couple to the three-dimensional gluon 
fields via the second term in Ti. shifting the longitudinal momenta of the trans- 
verse gluon fields (Eq. (|39|)). We neglect these couplings and consider the pure 
zero mode Hamiltonian 



h = I d 2 



x 



±.p 3 J(x ± )p 3 _(x ± ) + J^(V x a 3 _(x ± )f 



(40) 



We recognize 'electric' and 'magnetic' contributions in h, the zero mode Hamil- 
tonian - the first and second term, respectively. The light front variables mix 
the ordinary spatial and time variables so the labeling above is to be understood 
in analogy with the equal time Hamiltonian. Even at this level of approxima- 
tion, this zero mode Hamiltonian differs from the corresponding one in QED. 
The reason is the hermiticity defect of the canonical momentum: p_ ^ p_ . 
We now omit the color index and introduce dimensionless variables 

<p(x±) = ^a-(x x ) . (41) 

In the Schrodinger representation we then obtain 



d x± 



(42) 

where J(<p) is the Jacobian and equals the Haar measure of SU(2) 

JMx±))=sm 2 (v(x ± )). (43) 

The Jacobian is connected to the hermiticity defect of p- ; which stems from 
the gauge fixing procedure taking into account the curvilinear coordinates. 
The measure also appears in the integration volume element for calculating 
matrix elements. As in earlier approaches E9, ip will be treated as a compact 
variable, < tp < ir. 
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At this stage it is necessary to appeal to the physics of the infinite mo- 
mentum frame to factorize the reduced true energy h TC d and the Lorentz boost 
factor = i , since essentially h is a light front energy, and it is well known 
how it behaves under a Lorentz transformation. Thus we rewrite ft. rcc i = 277/1 . 
Since the integral over transverse coordinates can contain arbitrarily small 
wavelengths, we regularize h rc d by introducing a lattice to evaluate the trans- 
verse integral. The lattice vector b numbers the lattice sites, and e*i and e-i 
are the two unit vectors on the two-dimensional lattice. In order to have 
standard commutation relations on the lattice the derivative on the lattice be- 
comes — Sip ^ x± ^ a 2 . We further explicitly pull out the dependence on the 

lattice cutoff by defining a new reduced Hamiltonian /i rcc i = ait, d = 2rjah and 
substituting 77 = -j^aM, where M is a typical hadronic masscil. This yields 

ied = (44) 
b b 

with the electric contribution 



and the magnetic term 



hf = ^J2(n-n + f- (46) 



Since the effective coupling constant, 



2 g 2 LM 

5cff = 17T' (47) 



contains the large factor LM, the product of lattice size in the longitudinal 
direction and the hadron mass, a strong coupling approach seems to be a 
good starting point. Note that we avoid introducing 'radial wave functions' or 
effective potentials as others have done E3o. For each lattice site b, h% (the 



kinetic energy) has the Gegenbauer polynomials C n -(ip-g) for eigenfunctions: 
"Vpn~M) = Scff^K + 2)C nf fo> 5 ), (48) 

with 

b V 7T l_ sinyg J 
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and 

/ J(<p)C n (ip)C m {(p)dip = 5 nm . (50) 
Jo 

The strong coupling wave functions of the full transverse lattice are product 
states characterized by a set of quantum numbers {n} = {fir}, 

* W M=ii c « f K-)- ( 5i ) 



These functions form a complete and orthonormal basis for the zero mode 
sector. They satisfy the energy eigenvalue equation 

Y,h^{n}^)=9 2 cS J2 n b( n b + 2 )^{n}( l P)- (52) 



The ground state in this limit corresponds to all = - a constant wave 
function 



*{o}(^)=nvf ; 



(53) 



and the ground state energy is zero 

E Q = 0. (54) 

The first excited energy level is iV^-fold degenerate - an excitation at a single 
lattice point 

^^JI^M^JI. (55) 

b'^b 

In strong coupling this level is separated by a large amount from the ground 
state energy 

Ei = 3g 2 cS . (56) 

So far our results are equivalent to those reported by others El to within re- 
definitions of wave functions and integration measures. Weak coupling varia- 
tional solutions for the full SU(2) lattice Hamiltonian have also-been givenl23. 
Furthermore, studies in (l+l)-dimensional Yang- Mills theory c3 give formal 
extensions to construct wave functions for SU(N) gauge theories. 

The magnetic term of the Hamiltonian couples nearest neighbor lattice 
points. In the strong coupling limit its contribution may be obtained pertur- 
batively (as it has the coefficient l/g^g) by evaluating it with the basis function 
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of the ground state. The result of this is 

<*{0}lX>fl* { } ) = ^ (y - l) • (2iVi) . (57) 

Since this energy is proportional to N]_ = (L/a) 2 and g 2 s grows linearly with 
L, this part of the zero mode dynamics represents a negligible surface effect 
for the three-dimensional system in the strong coupling approximation. 

Next, we discuss the weak coupling limit g 2 — > 0. In this case we can 
simplify the kinetic term of the Hamiltonian by defining new variables a^: 

a,= f g , (58) 
with 8k 2 = \J~2LM . After expanding in g, the reduced Hamiltonian becomes 

^ = E {- + ^) + Ek - a ^ r } • (59) 

Going over to Fourier momentum representation, 



E e ^%> ( 6 °) 



with ki = 2nrii/N±a and = 0, ±1, ±2, ... we have 

»- - E {- (|| + ^4) + 4 X> 2 IV-t} ■ <«> 

The eigensolutions tAa: of h le d in the weak coupling approximation are decou- 
pled harmonic oscillators for each k, with frequencies 

<4 = 4 £ sin2 y- (e 2 ) 

e 

Because of the 'radial Laplacian' it looks as if the eigenfunctions would have 
to vanish at the origin to be normalizablc. However, as in the Schrodinger 
equation in three dimensions, the Jacobian J allows a constant wave function 
at the origin. Consequently, the eigenvalue of ipx is given by the sum over the 
modes: 



n* = E* 



4 E sin2 (¥)' ( 63 ) 
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which gives in the N± — > oo limit spin waves with = \Jk\ + fcf . 

In the weak coupling limit the zero mode Hamiltonian supports solutions 
similar to QED. The strong coupling limit, however, yields different results: 
'gluonic' excitations are suppressed because of large energy gaps. This is due 
to the Jacobian, which can be traced back to non-abelian self interactions in 
the original Lagrangian. 



4-3 Effective Hamiltonian for Two-Site Truncation 

We will now invoke an effective interaction approach embedded in a cluster 
expansiontil, starting with the simplest, two-site cluster, in which either site (or 
both) can be excited to high energy states. We will obtain the solution for the 
low-lying spectra of the system approximated as a low density of excitable two- 
site clusters over the entire range of coupling. This method can be envisaged 
as the starting point of more ambitious renormalization group techniquesli-i ' . 

We work in the representation of the strong coupling solution of /i2 and 
divide the two-site subspace into a P and Q space, such that P + Q = 1 with 

P = {|0,0)}, 

Q = {|n,m); n,m ^0,0}, (64) 

where n, m represent the indices of the Gegenbauer polynomials. Then the 
two-site energy E% is given by the non-perturbative solution of the Hamiltonian 
in Eq. (Q), truncated to two lattice sites and labelled h 2 . Explicitly, this 
Hamiltonian is: 

ha = h c + h m , (65) 

with 

Scff | j g^ + j fi^ | ' ( ) 

and 

fi m = - W)\ (67) 

where the subscripts label the sites. 

i_ Within the effective Hamiltonian method, the two-site energy is given 
byN: 

E 2 = Ph 2 P + PhaQ — Qh 2 P . (68) 

E 2 - Qh 2 Q 

The self-consistent solutions of this equation provide the low-lying spectra in 
this method. The strong coupling basis states are eigenstates of h e : 

h c \n, m) = g 2 cS {n(n + 2) + m(m + 2)} \n, m) . (69) 
15 



Thus, the non-trivial matrix elements are those of h m , and are of the form 
(n,m\(ipi — Lp 2 ) 2 \n' ,m') , where the subscripts label the sites. The matrix 
element reduces to sums and products of one-site matrix elements, which are 
given as: 



{n\ip\n ) 
(n\(p\n') 



l( i 

7r \ (n+n' 







(n-n>)* 



for n = n 

for n + n' = odd 

for ri + n' = even, n =/= n' 

for n = n' 



l l 

(n-n') 2 (n+n'+2) 2 



| for n =/= n! . 



(70) 



In the strong coupling limit we obtain the result of perturbation theory in 

l/Seff 



E 2 = (0,0|n0,0) = -L {(0|^ 2 |0) - 2 «0|^|0)) 2 } = i- (lL - 1 



(71) 



In the weak coupling limit we can solve the Schrodinger equation for the 
two-site version of the earlier spin wave Hamiltonian, Eq. ([59]). We call the 
respective variables ctg — x and = y , then solve: 



(—- 


2 d \ 


(£.. 


\- 2 9 \ 


\dx 2 


x dx J 




y dy) 



+ (x- yf 



(72) 



Since this Hamiltonian is invariant under x <-* y, the eigenfunctions ^2{x,y) 
can be chosen to be symmetric under the interchange of x and y (^2s(x, y)), or 
antisymmetric (Jl l 2a{x,y)). The first symmetric excited state becomes degen- 
erate with the ground state of the original problem in the weak coupling limit, 
and the first antisymmetric state has a greater energy than the first symmetric 
state. 

As usual one factorizes the wave function ^(^ y) — —&2(x, y) , resulting 
in the Schrodinger equation 



if 



d 2 



(73) 



The center-of-mass motion is then separated: $2(2;, y) = e Xi{ r ) > with 
R = (x + y)/2 and r — x — y. The Hamiltonian corresponding to the relative 



1G 



motion (r) is a simple radial harmonic oscillator: 

(^ed), = -2^+r 2 . (74) 

The lowest states of the symmetric and antisymmetric 'towers' are solutions to 
this Hamiltonian. The energies of these states can be read directly from Eq. 
([74]); E 2s = V2, and E 2a = 3\/2, respectively, giving a energy gap between the 
states of 2^2. 

Thus, the results for the energy gaps of the low-lying states in the weak 
coupling limit are 

E 2s — Aground = , 

E 2a — Aground = 2v2 . (75) 

We now proceed to calculate the low-lying spectra via the effective Hamil- 
tonian method. In the numerical calculations, we truncate the Q space at a 
certain two-site energy, calculate E 2 , then increase the energy until we reach 
convergence for each choice of coupling constant g e ff- The typical number of 
two-site states required for the Q space at convergence was 300. ._. 

The numerical solution of Eq. ( |68| ) for E 2 has now been publishedtZI. In the 
strong coupling limit (large g^g) the large gaps in energy are evident, and the 
numbers agree with the unperturbed energy of the states, given in Eq. (|6^) . In 
this same limit, the slope of the ground state energy as a function of the inverse 



square coupling agrees with the analytic calculation of Eq. (71). In the weak 
coupling limit, the results for the gap energies were Richardson extrapolated 
for the l/<7gfi ~~ > 00 mlu t- This extrapolation matched the analytical results of 
Eq. (|75| ) to five significant figures. Thus, we have obtained two-site solutions 
in the strong coupling basis for the entire range of coupling which agree with 
analytic results in the both the weak and strong coupling limits. 

It is remarkable that we succeeded with a one-dimensional strong coupling 
basis throughout the range of coupling strengths. Results for the spectra of 
iV-sites are straightforwardly obtained and should be valid as long as the 
number of excited two-site clusters is small compared with N/2. This is the 
'low density' approximation. 

5 Variational Tamm-Dancoff Application 

Over the past several years, we have also investigated the use of relativistic 
wave equation formulations of quantum field theory in the equal time quan- 
tization scheme. Here, we would like to briefly summarize our approach as 
viewed from the effective Hamiltonian framework. E§Ea 
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We begin by selecting the Coulomb gauge Hamiltonian of QCD and we 
choose to work in a finite momentum space domain. We eliminate all interac- 
tions except the fermion-gluon vertex and the triple gluon coupling terms. We 
fix the scale by giving the quarks a constituent quark mass. 

For our unperturbed Hamiltonian we select the kinetic energy operators 
for the quarks and the gluons. In our earlier workcJ we chose a cubic B-spline 
basis expansion for our quark and gluon amplitudes. More recently,c£l we select 
a harmonic oscillator basis with a variable oscillator parameter, resulting in a 
non-orthogonal basis. 

For our P-space, we select a small finite set of quark-antiquark states for 
mesons and 3-quark states for baryons. 

Our procedure is to optimize our P-space using a variational approach and 
to determine an effective interaction for that basis. The effective interaction 
emerges as a ladder series as given by the G-matrix defined earlier. However, 
there are a number of simplifying assumptions needed to evaluate this effective 
interaction. 

Finally we solve the resulting Eq.(||) which is the Bethe-Salpeter equation 
with the kernel defined by the effective interaction we have derived. For fixed 
constituent quark masses and a fixed strong coupling constant., the solutions 
of this equation yield spectra in reasonable accord with data. 

The most remarkable feature of this approach is found in the shape of the 
effective interaction which emerges. At short distances it closely approximates 
the one gluon exchange effective potential and, at large distances, it closely 
approximates linear confinement. 

The linear confining feature is traced to the non-linear role of the triple- 
gluon coupling term and its effect in our variational treatment. The combina- 
tion of the choice of Coulomb gauge and the variational approach appear to be 
ultimately responsible for this result. We believe that this feature will persist 
as we make our approach more general. 



6 Conclusions 

The well established framework of effective Hamiltonians for quantum many- 
body theory with strong interactions appears adaptable to non-perturbative 
applications in quantum field theory. Our initial results on simple renormaliza- 
tion problems are encouraging and the initial applications to (simplified) QCD 
problems also inspire us to proceed. However, we realize that considerable 
effort is required for major additional advances to be achieved. 
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